An alternative approach to exact wave functions for time-dependent 
coupled oscillator model of a charge in variable magnetic field 

Salah Menouar 1 , Mustapha Maamache 1 and Jeong Ryeol Choi 2 

1 Laboratoire de Physique Quantique et Systemes Dynamiques, 
Departement de Physique, Faculte des Sciences, 
Universite Ferhat Abbas de Setif, Setif 19000, Algeria 
(E-mail: menouar_salah@yahoo.fr) 

2 School of Electrical Engineering and Computer Science, 
Kyungpook National University, 1370 Sankguk-dong, 
Buk-gu, Daegu 702-701, Republic of Korea 
(E-mail: choiardor@hanmail.net) 

Abstract 

The quantum states of time-dependent coupled oscillator model for charged particles 
subjected to variable magnetic field are investigated using the invariant operator methods. To 
do this, we have taken advantage of an alternative method, so-called unitary transformation 
approach, available in the framework of quantum mechanics, as well as a generalized canonical 
transformation method in the classical regime. The transformed quantum Hamiltonian is 
obtained using suitable unitary operators and is represented in terms of two independent 
harmonic oscillators which have the same frequencies as that of the classically transformed 
one. Starting from the wave functions in the transformed system, we have derived the full 
wave functions in the original system with the help of the unitary operators. One can easily 
take a complete description of how the charged particle behaves under the given Hamiltonian 
by taking advantage of these analytical wave functions. 

Key Words: Time- dependent systems; coupled oscillator; unitary transformation 
Schrddinger equation. 
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1. INTRODUCTION 

Since the pioneering works of Lewis[l-3|, the investigation of the quantum behavior for time- 
dependent harmonic oscillator has attracted considerable interest in the literature because it 
offers exactly soluble models for dynamical systems in different areas of physics. There are 



a a, 



diverse kind of time-dependent harmonic oscillators such as Caldirola-Kanai oscillator 
parametric oscillator jf]], and harmonic oscillator with a strongly pulsating mass The two- 
dimensional time-dependent harmonic oscillator also became an hot topic in both classical 
and quantum mechanics. The higher dimensional harmonic oscillator has played an important 
role in, for examples, nuclear shell structure and models of quark confinement. According 
to the progress of research for these systems, a relevant question has been naturally raised: 
what would happens if two-dimensional harmonic oscillator is constructed by the couplin g of 
external additive potentials? The investigation for this issue was initiated by Kim et al.j8-3] 
about thirty years ago. They considered two harmonic oscillators that are mutually-coupled 
so that the resulting potential becomes V(xi,x 2 ) = \ [c\x\ + c 2 x\ + 03X1X2). They determined 
the corresponding density matrix rigorously in order to establish the Wigner function and 
some of other useful functions in physics. There are plenty of physical systems and mode 



described by coupled harmonic oscillators, such as the Lee model in quantum field theory 14], 
the Bogoliubov transformation model of superconductivity 11511 . two-mode squeezed light [161 ]. 
the covariant harmonic oscillator model for theparton picture jl7|. and some models in molecular 



physics 



181 ] . One of the main focuses of research fulfilled in connection with time-dependent 



coupled oscillators is some specific problems of time-dependent coupled electric circuit whose 



closed-form solutions are now well known 



19 



221] . And further, the author of Ref. 1 231 have 



investigated the propagator for a certain class of time-dependent coupled and driven harmonic 
oscillators with time- varying frequencies and masses using path integral methods. 

Above all, a charged particle in a strong uniform magnetic field is a typical example of natural 

n 

non-commutative system |24|, which provides a good starting point when we are going to discuss 
the quantum Hall effect [25|]. The external magnetic field is indeed an important factor that 
affects to the motion of a charged particle in various physical systems. When a time-dependent 
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magnetic field is exerted on an electron, it is impossible to reduce the system to be a one- 
dimensional problem. Instead, it can be modeled by two-dimensional time-dependent harmonic 
oscillator due to the existence of variable magnetic field B(t). Theoretical and experimental 
researches have been carried out extensively on the quantum properties of this system in the 
past several decades due to its importance not only in condensed matter physics but also 



in plasma physics 



26- 



Especially, the study of charged particle motion driven by external 
magnetic field is crucial in investigating magnetic confinement devices for fusion plasmas (whose 
subtopics are tokamaks, mirror machines, bumpy tori, and stellarators) and for space and 
astrophysical plasmas (whose subtopics are magnetospheric plasmas of earth and other planets 
like pulsar) 36]. The wave functions of a free electron with time-dependent effective mass, in 
;he presence of a variable magnetic field, are derived in both the Landau and symmetric gauges 



37j |. The propagators of a charged particle subjected to a time-dependent magnetic field, which 



propagate the wave functions in the spacetime, are derived using the linear and the quadratic 
invariants [38j] . 

An interesting problem that would be worth to be dealt with is the charged particle sys- 
tem that is described by the Hamiltonian which involves the static coupling xy and dynamic 
coupling term p x p y under the presence of magnetic field. This system may exhibit novel fea- 
tures owing to the existence of the coupling terms and can be regarded as the generalization 
of the Hamiltonian model given in Refs. 31] and 39) . Our intention in this paper is to cal- 
culate the exact wave function for time-dependent coupled oscillators in a variable magnetic 
field within the framework of the invariant methods. The calculation is based on the use of 
the generalized time-dependent canonical transformations and an alternative time-dependent 
unitary transformation. 

The present paper is organized as follows. Our problem is formulated in section 2 through 
a general time-dependent Hamiltonian describing the motion of a complicate charged particle 
system. Some remarks necessary in dealing with our task will also be presented. In section 3, 
we show how to simplify the problem associated with the complicate Hamiltonian of our system 
using the canonical transformation method. As an alternative approach, unitary transformation 
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is also applied, in section 4, in order to transform our complicate Hamiltonian to that of a more 
simplified harmonic oscillator. The quantum solution of the system will be investigated in 
section 5 using the invariant methods on the basis of the results obtained in section 4. The 
concluding remarks are given in the last section. 



2. FORMULATING THE PROBLEM 



Let us formulate our problem by introducing a generalized Hamiltonian describing the motion 
of a charged particle that have time-dependent effective mass m(t) in the presence of a variable 
magnetic field. The effective mass of charged particles, such as electrons or holes in any system, 
may modified through their interaction with surroundings or various excitations like energy 40] , 
stress 41], temperature 42], and pressure [431]. It is therefore natural to think that the effective 
mass varies with time according to the change of the environments. Moreover, if we vary the 
external magnetic field randomly in the heterojunctions and solid solutions, the effective mass 
of an electron also varies in a random fashion in response to the fluctuation of the composition 
in the system [441]. 

We consider the electromagnetic potential in the symmetric gauge such that 
A ( — f^ y, ^p-x, 0) where x and y are the position operators. Then, for the dynamical sys- 
tem of our interest, the Hamiltonian has the form 

H(x, y, t) = (n 2 x + ttJ) + l -m{t)w 2 {t) (x 2 + y 2 ) + a{t)xy + b(t) PxPy , (1) 

where ir x = p x — eB ^ y and ir y = p y + eB ^ x while p x and p y are conjugate momentum operators 
that are given by p x = —iTid/dx and p y = —ihd/dy. zn(t) is the oscillating frequency that is an 
arbitrary function of time. To generalize the problem, we suppose that the other parameters, 
m(t), a(t) and b(t), are also arbitrary time functions. 

In terms of p x and p y , the Hamiltonian (1) can be rewritten as 

H(x, y, t) = (pi + p 2 y ) + l -m{t)uj 2 {t) (x 2 + y 2 ) + a(t)xy + b(t) PxPy + ^p-L z . (2) 

Here, u(t) is a modulation frequency which takes the form co 2 (t) = w 2 (t) + where w c (t) = 

4 



is the Larmor frequency, and L z = xp y — p x y = —ih (x-^ — y-§^j is the canonical angular 
momentum in the axial z direction. The study of quantum features for this system may be a 
fascinating task in its own right both from a physical and a mathematical point of view. We will 
show how to convert this Hamiltonian to a simple form in the following two sections by means 
of the canonical transformation and unitary transformation, respectively. These procedures 
may enables us to derive the quantum solutions of the system. 

3. CANONICAL TRANSFORMATION 

The method of time- dependent canonical transformation is in fact very powerful in inves- 
tigating the mechanical behavior of dynamical systems. We can convert a given Hamiltonian 
into any desired one through this method. In order to cast the Hamiltonian of our problem into 
a more soluble form, we take the advantage of the time-dependent canonical transformation 
(x, y, p x , p y ) (X, Y, P x , P Y ) defined as. 



x\ _ I cos 00) sin 00) \ fX 
v) \-sm(j){t) cos ) V F 

p \ = I cos0(t) sin 00) \ (P x 
Py) \ - sin <f>(t) cos <p{t) ) \ P y 



(3) 



(4) 



where 

m = -\J ^c(t)dt. (5) 
From the fundamentals of classical mechanics, we have |fl 

Vx = —F 1 {x.y,P x ,P Y ,t) ,X = —-F 1 (x.y,P x ,Py,t), (6) 
ox oP x 

Py = §~ F i fay, Px, Py, t) , Y = -^-F 1 (x.y, P x , P Y , t) , (7) 



and 



(P X X + P Y Y - H(X, Y, t) = Px x + p y y - H(x, y, t) + (8) 
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where i*\ is the generating functions responsible for the transformation. Through these rela- 
tions, Fi is easily found to be 

Fi (x.y, Px, Py, t) = (Px cos + Py sin 0) x + (— Px sin + Py cos 0) y, (9) 

= -0(i)L 2 = ^L 2 . (10) 

In terms of the new conjugate variables (X, Y, Px, Py, ), the Hamiltonian (2) becomes 

H(X,Y,t) = ^^ + ^^ + ^-(^(0^ + ^ + (^K(^ 2 

+a 1 (t)XY + b 1 (t)P x Py, (11) 

26(t)sin0cos0, (12) 

26(t) sin cos 0, (13) 



where 



m-(t) m(t) 
1 1 



m + (t) m(t) 

1 /2 

f m(t)uj (t) — 2a(t) sin0cos0\ A , 

= ( *,.(*) J ' (14) 

= ^ m(t)^(t) + 2a(t)sin^cos^ 1/2 

ai (i) = a(f) (cos 2 - sin 2 0) , (16) 

6i(t) = b{t) (cos 2 - sin 2 0) . (17) 

To eliminate the dynamical term PxPy, we take the second canonical transformation by re- 
casting the canonical variables (X, Y, Px, Py) in terms of new variables (gi, qi, Pi, P2)'- 

X = (!^y\ " 4 (W) , Y = (^m) " 4 + , (18, 



X - 7TT , r Y 



The canonical transformation based on Eqs. (18) and (19) enables us to transform H(X,Y,t) 
into H(qi,q 2 ,t). Thus, straightforwardly, we have 

F(gi,g2,t) = 2^t/ 1 + 1 hhj/ 2 + l 7711 ^^ 1 + ^ m 2(t)w 2 2 (t)g 2 2 + c(t) gi g 2 , (20) 



where 



1 



1 



m^t) y/m+(t)m-(t) 

1 _ 1 
m 2 (t) y/m+(t)m-(t) 



-h(t), 



w 2 (t) 



Vm + (t)m_(t) [(<«£(*)+ <4(t))] -ai(t ) 
mi(*) 

|y/m + (t)m_(t) [(a; 2 (t)+a; 2 (t))] + ai (t) 
m 2 (t) 

it) = ±y/m + (t)m_(t) (u 2 _(t) - ^(f)) . 



1/2 



1/2 



(21) 

(22) 

(23) 

(24) 
(25) 



Now, to remove the static coupling term qiq2, we take anot 



rer canonical transformation by 



introducing the variables (Q i7 Pi) where i = 1,2, such that [20 



21 



23] 



1 9(t) 
, COS -rr 

yWW 2 



v/mi(t) 2 



1 • 0(f) 1 

— i — sin — - — 

A/ m 2(*) 2 a/^W 



\Jmi{t) cos -4^ ■s/mi(t) 



sin 



0(f) 



cos 



Pi 



0(f) 



/ 



Qi 

Q 2 




(26) 



(27) 



-\/ m 2(i) sin ^ ^m 2 {t) cos / 
where 0(f) is an arbitrary phase which shall be appropriately determined afterwards. Equations 
(26) and (27) do not always represent canonical transformation 45j between the variables (qi,Pi) 
and variables (Qi,Pi). If (Qi,Pi) are canonical coordinates, there should exist a new Hamilto- 
nian which is determined only by the Hamiltonian given in (20) and the linear transformation 
of (26) and (27). 

The variables (qi,Pi) and (Qi,Pi) in the two representations must satisfy the following 



relation 



BF 

J2 p iQi - H Q = l>fc - H « + -Qr> ( 28 ) 

i=l i=l 

provided that the transformation is canonical, where F is called a generating function, possibly 



a time-dependent function in phase space. From the equations known in classical mechanics 
Pi = ^-F(q 1 ,q 2 ,P 1 ,P 2 ,t) ,Q i = -£ T F(q 1 ,q 2 ,P 1 ,P 2 ,t) 2 = 1,2, (29) 



the generating function responsible for the transformation is found to be 



F{q 1 ,q 2 ,P 1 ,P 2 ,t) = y/m^t) ^ cos ^ + P 2 sin Ql 

+\/m 2 {t) (-Pi sin d -^- + P 2 cos 



0(t) 



<l2 



+^[-^rh 1 (t)Q 2 i-^rh 2 (t)Q 2 2)- (30) 

In terms of the new conjugate variables (Qj,Pj) the Hamiltonian of the system can be 
rewritten as 

H Q (Q 1 ,Q 2 ,t) = ±(p? + P?) + ±nl(t)Ql + ±nl(t)Ql 

+ 6 ~Y- [P1Q2 - P2Q1] + 5(t)QiQ2, (31) 

where the time- dependent coefficients Qi(t),Q 2 (t) and S(t) are given by 

\ 1/2 

w-w^m+w^m-jp*™) , (32) 



.-.2^ • 2*(*) , ~2^ 2^(t) , c(*)sin0(i) 



1/2 



n 2 (0 = ( sm 2 -y + cos 2 -y + w v ; | (33) 

6(t) = \ (d>?(t) - ^ 2 (t)) sin(9(t) + f ] ™ e ® (34) 

with 

^ {t) = " l(t) + 4^^)- 2 ^)J' (35) 

-lit) =c, 2 (t) + if4|-2^V (36) 
4 V^2(t) m2(t) / 

If we choose 6*(t) = Const., the term involving P\Q 2 and P 2 Q\ in Eq. (31) disappears and, 
consequently, we have 

H Q (Q 1 ,Q 2 ,t) = l - (P 2 + P 2 ) + \n\(t)Q\ + l -tf 2 {t)Q 2 2 + (J(t)Qig 2 . (37) 
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It is notable that, with the above canonical transformation, the coupling S(t) is a function 
of the parameters of the original system. Evidently, the separation of variables in equation (37) 
requires that S(t) = , i.e. 



and consequently 



c(t) = ^ (w|(t) - ^ 2 (t)) v/mi(t)m 2 (t) tan0, (38) 

ta ^VW(«). (39) 

(w 2 2 (f ) - (t) ) 

Under this condition, the Hamiltonian of Eq. (37) reduces to 

H Q (Q u Q 2 ,t) = \ {Pi + Pi) + \n\{t)Q\ + \nl{t)Ql (40) 

This is the sum of two individual Hamiltonians corresponding to the harmonic oscillators having 
the time-dependent frequencies VL\{t) and fl2(t), respectively, and having masses that are equal 
to unity. 



4. UNITARY TRANSFORMATIONS 



The unitary transformations in quantum mechanics is analogous to the canonical transforma- 
tions in classical mechanics. In this section, the relationship between the two transformations 
will be demonstrated and we confirm how to obtain the quantum-mechanical Hamiltonian from 
the classical one. With the consideration of quantum physics, we replace the canonical variables 
(x,y) by quantum operators (x,y), so that the corresponding Hamiltonian is given by 

H(x, y, t) = (pi + pi) + l -m(t)uj 2 {t) (x 2 + f) + a(t)xy + b(t)p x p y + ^-L z . (41) 

In fact, it is not difficult to show the commutation relations [L z , x 2 +y 2 ] = and [L z , Pl+p 2 ] = 
0. We can also check the non-commutability of L z with some other variables: [L z , p x p y ] ^ 0, 
and [L z , xy] ^ and consequently [L z , H] ^ 0. Considering this fact, we are unable to 
decompose the Schrodinger equation 

ih^(x,y,t) = H(x,y,t)V(x,y,t), (42) 
9 



when we would like to simplify it, because the angular momentum operator L z and the Hamil- 
tonian H(x, y, t) does not have the same eigenstates. To overcome this difficult situation, we 
transform the Hamiltonian (41) to a simple form by means of appropriate unitary operators. 
In the first step, we perform the unitary transformation 



*(x,y,t) = U(t)ij>((x,y,t), 



(43) 



where U (t) is a unitary operator of the form 



U(t) = exp 



' 2h 



I 



zu c (t)dt 



(44) 



Under this transformation, the Schrodinger equation of original systems (42) is mapped into 



d 



(45) 



where the new Hamiltonian Hi(x,y,t) has the form 



Hi(x,y,t) = jjj fil + 2m + (tf 2y + ^ m -(*) w -(*) £2 + ^ m +( t ) UJ l( t )f + ai(t)xy + bi(t)p x p y . 

(46) 

Note that the term involving L z disappeared in the above equation, This means that the mag- 
netic field is removed when it is viewed from an appropriate rotating frame (p(t) = — | J w c {t)dt. 

To simplify the Hamiltonian (46), we use two-step unitary transformation approach. As a 
first step, we take the following unitary transformation 



where 

and Ai(t) and A 2 (t) are given by 
A x (t) = exp 



$(x,y,t) = A(t)ip(x,y,t), 



A(i)=A 1 (t)A 2 (i), 



-( Px x + x Px )ln{— - 



exp 



, , m + (t)\ 



1/4' 



k 2 (t) = exp 



l 7T 

hi 



(xp y - yp x ) 



(47) 



(48) 



(49) 



(50) 
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Then, we can transform the Hamiltonian (46) using the formula 

H 2 (x, y, t) = k~ l {t)Hi(x, y, i)A(f) - ihk~ l (t)^-k{t). 

After some algebra, we get 

1 



(51) 



H 2 (x,y,t) 



:P X + 



Py + ^MtMiW + -m 2 (t)oo 2 (t)y 2 + c(t)xy. (52) 



2 mi (ty x 2m 2 (ty y 2 ■" ir '"' lvv " ' 2' 

In the next transformation we will eliminate the coupled static terms xy. To do this we consider 
the unitary transformation 

<p(x,y,t) = V(t) X (x,y,t). (53) 
Here, V(t) is a time-dependent unitary operator of the form 



V(t) = V 1 (t)V 2 (t)V 3 (t), 



where 



V\{t) = exp 



2h 



(p x x + xp x ) In y/mi(t) 



exp 



2h 



(PyV + VPy) In y/m 2 (t) 



V 2 (t) = exp 



-y^Py-VPx) 



V 3 (t) = exp -— (m^x 2 + m 2 (t)f) . 



(54) 

(55) 
(56) 
(57) 



Substituting equation (54) in equation (53), we can obtain a transformed Hamiltonian that is 
merely the coupling of two harmonic oscillators having frequencies fli(t) and fl 2 (t) and unit 
masses: 



H 3 (x,y,t) = V-\t)H 2 (x,y,t)V(t)-ihV-\t)^V(t) 

= \{pi+Pi) + \m- 2 +\mf. 



(58) 



At this stage, one can easily confirm that the relation given in equation (40) is correct, since 
it is consistent with equation (58). From unitary operators (44), (49), (50), (55) and (56), we 
can confirm that Ai(t) and Vi(t) are the squeeze operators whereas U(t), A 2 (t) and V 2 (t) are 
the rotation operators with the angles 0(t), | and |, respectively. 
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Though the original Hamiltonian (41) involves the static coupling term xy and the dynamic 
coupling term p x p y , the transformed Hamiltonian (58) does not have such terms. Hence we 
can easily handle equation (58). In the following section, we establish the quantum solution 
(wave function) in the transformed system. And then, we will take the advantage of the unitary 
transformation starting form this wave function using the same unitary operators introduced 
in this section in order to derive the full wave functions in the original system. 



5. QUANTUM SOLUTIONS 

The problem of the harmonic oscillator with time- dependent mass and frequency can be 
transformed to that of the harmonic oscillator via the associated invariant [3j. It is easy to 
verify from Liouville-Von Neumann equation 



dI 91 1 rf frl n 



(59) 



that the invariant associated with the transformed Hamiltonian of two-dimensional harmonic 
oscillator is given by 



I(x,y,t) = I(x,t)+I(y,t) 

y^j + {pi% - p\x)' 



1 

+ 2 



— ) +(p2y-p2y)' 

P2, 



where pi(t) and pzif) are c-number quantity satisfying the auxiliary equations 

P \ + n 2 1 {t) Pl = l/pl, 

p 2 + Ql(t)p 2 = 1/pl 



(60) 



(61) 
(62) 



In order to make the invariant hermitian, P = I, we choose only the real solution of (61) and 
(62). Then, we can find the eigenfunctions £, nin2 (x,y,t) of I(x,y,t), which are a complete 
orthonormal set corresponding to the time-independent eigenvalues A„ iri2 , from the eigenvalue 
equation 

I(x,y,t)£, nim (x,y,t) = \ nin2 £ niri2 (x,y,t), (63) 
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where 



Kn 2 = % (ni + 2 J + h ("2 + 2 J • ( 64 ) 
By evaluating equation (63) with the use of equation (60), we have the eigenstates in the form 



£nin 2 (*£) 2/) 



1 



7r^ni!n 2 !2 n i+ n 2pip 2 
« / pi , i 



1/2 



x 



x 



ni U 1/2 



pi 



V 



h 1 ' 2 



pi 



x exp 











\P2 


pi) . 



_2h Vpi p\, 

where H ni and H n , 2 are the usual Hermite polynomial of order ri\ and n 2 , respectively. 
The solution of the Schrodinger equation 

d 

ih—x ni n 2 (x,y,t) = H 3 ( y x,y,t)xn 1 n 2 (x,y,t), 
can be written in the form 

Xnin 2 Vi t) — e nin2 ^C, nin2 {x,y,t), 

where the phase function a nin2 (t) satisfy the equation 
d 1 d 

-fo a nm 2 (t) = (Cmna J/> *) I ^ ~ #3(z-J/, |£mn 2 J/, *)) ■ 



(65) 



(66) 



(67) 



(68) 



According to equations (67) and (68), the solutions Xnm 2 (%, V, t) of the transformed Schrodinger 
equation (66) are given by 



V (t 11 f} = <s* an l n 2(*) 

x exp 



1 



7ihn 1 \n 2 \2 ni+n2 p 1 p 2 



1/2 



X Hr, 



X 



111 Vft 1/2 



pi 



H„ 



^ 1/2 P2 



Pi 



2ft U + P?r + 2ft U + Pi, ,r 



P2 



where the phase functions takes the form 



1 dt' ( 1 

7T " n 2 + o 



dt' 



(69) 



(70) 



V Jo pW) 

The relation between the wave functions, ty nin2 (x, y, t) in original system described by the 
Hamiltonian (41) and the wave functions Xn in2 (x,y,t) in the transformed system is 



* nina (x,y,t) = U(t)A(t)V(t) X (x,y,t) 

= f>(t)A 1 (t)A 2 (t)V 1 (t)V 2 (t)V 3 (t)Xn 1 n 2 (x,y,t). 



(71) 
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Using equations (69), (70) and (71), we derive the full wave functions in the form 



1/2 



1/2 



(mi(t)m 2 {t)) 

(rji cos + 772 sin 0) x + (—771 sin (ft + i] 2 cos 0) y 
(/iii cos + /i2 sin 0) x + (— //1 sin + p 2 cos 0) 7/ 



xH, 



in 



n^ P2 



xexp 2h 



x exp 1 



X 1/2 



m 



m 



1/2 



1/2 



j /1 cos 2 + ( — — ) f 2 sin 2 + ^ sin 20 



/icos^0+ ^ / 2 sin^0-^sin20 

m_i_ / \ 77i- / 2 



1/2 



fs . 



— / 2 - — /1 ) sin 20 + / 3 cos 20 

77l_ / V 777_ / 



1 

^1 + 77 



where the time-dependent coefficients / 2 (0> /3(0> 7(0 , /3(0> 771(0, rj 2 (t) : pi(t), and /x 2 

are given as follows 



* dt' 



/i(0 = + ^m 2 ^j cos 2 0/2 + Qm 2 + |mi ) sin 2 0/2 

J r^/mim 2 (/3 — 7) sin0/2 cos0/2, 



(■ 



/ 2 (0 = (|mi + 7^2) cos 2 0/2 + Qm 2 + ^ ) sin 2 0/2 
—y / m 1 m 2 (/3 — 7) sin 0/2 cos 0/2, 



/ 3 (t) = (-7777! + /3m 2 ) cos 2 0/2 + ( 7 m 2 - /5mi) sin 2 0/2, 

1 J_ _ ^i(0 

2 

vPl Pi 



7(0 



£(0 
\ 1/4 



P2 J_ _ 

P2 pi 



2 

m 2 (t) 



?7i(0 = (V™i/2cos0/2 - v/m 2 /2 sin0/2) , 



c 


c 

(■ 
c 
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/ s 1/4 

V2<t)=l?±) (-VW/2 cos 9/2 - ^m 2 /2 sin 0/2J , (79) 

1/4 



(V™i/2sm#/2 + y/m 2 /2cos6/2 s } , (80) 

/ \ 1/4 

/i 2 (t)=(^ ± J (- y /m 1 /2sm6/2+ y/m2/2cos6/2) . (81) 

The final solutions given in (72) are somewhat complicate, but they are very useful when 
predicting the evolution of the probability distribution of the system. A considerable physical 
significance of such analytical solutions is that their application in physical system is very 
flexible, even when the internal and external situations of the system vary from time to time. 
Numerical solutions obtained from, for example, the FDTD (finite difference time domain) 46] 
method are however inconvenient as inputs to further analyses, especially when the parameters 
of the charged particle system vary with time like in this case. One can easily take a complete 
description of how the charged particle behaves under the given Hamiltonian, by means of this 
analytical wave function. 



6. CONCLUSION 



Though the motion of charged particles in magnetic fields is a fascinating problem in both 
quantum and classical view, most of the relevant research is concentrated on static problem that 
can be described by time- independent Hamiltonian. In this paper, this problem is generalized 
to a more complicated case that is described in terms of time-dependent Hamiltonian, by 
supposing that the parameters such as the effective mass of the charged particle vary explicitly 
with time in the presence of variable magnetic field. We have presented an alternative treatment 
after reducing the problem related to the charged particle motion to that of the quantal time- 
dependent coupled oscillators. 

We approached the problem in two ways. In one, we used a time-dependent generalized 
canonical transformations which enabled us to transform the initial classical Hamiltonian (2) 
to a more simplified one associated with the two harmonic oscillators having time-dependent 
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frequencies fl\(t) and f^^)- We have taken, as an another way, an alternative approach based 
on unitary transformation that allowed us to transform the quantal Hamiltonian (41) to an 
equally simple one, but in the framework of quantum mechanics. To facilitate the derivation 
of quantum states, we introduced dynamical invariant operator I(x,y,t). The eigenvalues and 
eigenstates of the invariant operator are obtained using the Liouville-Von Neumann equation. 

We derived the exact wave functions of the system on the basis of the fact that they are the 
same as the eigenstates of the invariant operator, except for some time-dependent phase factors 
exp [ia ni „ 2 (i)]. The wave functions in the transformed system presented in equation (69) are 
relatively simple and expressed in terms of Hermite polynomial. However, as you can see from 
equation (72) with equations (73)-(81), the wave functions in the original system are somewhat 
complicate. These wave functions are represented in terms of pi and p 2 that are the solutions 
of some classical equation of motion given in Eqs. (61) and (62), respectively. If we let x C) i and 
x Cj 2 as the two independent classical solutions of the ^-component transformed Hamiltonian 
presented in (58), p\ can be written as p\ = [aix 2 cl + a 2 x Ci iX Ci2 + x^ 2 ] 1//2 (see Ref. 47] for 
rigorous mathematical proof). Of course, y-component two independent solutions y c> \ and y Cj2 
hold the similar relation with p 2 : p 2 = [biy 2 t + ^2?/ c ,i?/c,2 + ^y^ 1 ' 2 . This implies that we can 
take the complete knowledge for the behavior of the system within the scope that quantum 
mechanics admits when we exactly know the classical solutions of the transformed system, since 
we did not used any approxination in the development of our theory. 

Though we expressed the wave functions in terms of the classical solutions of transformed 
system, many authors represent the wave functions for time-dependent Hamiltonian system in 
terms of the classical solutions of original system. However, in this case, it is unclear that 
we can represent the wave functions using the classical solutions of original system since the 
classical solutions in the original system can not be decoupled into each coordinate component 
due to the existence of coupled terms xy and p x p y in equation (41). Even if it were possible to 
manage the problem using those of original system, the corresponding mathematical procedure 
to obtain the exact wave functions would be very difficult and the results would become much 
more complicate. The wave functions we obtained here can be used to evaluate not only the 
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quantum mechanical expectation values of various observables such as physical momentum and 
quantum energy but also probability densities and fluctuations of the canonical variables. 

Finally, we did not considered thermal effects in this work. The quantum behaviors of the 
charged particle motion with consideration of thermal effects may be a good topic as a next task. 
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